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Abstract 

Background: It is well known that mushrooms produce defense proteins and secondary metabolites against 
predators and competitors; however, less is known about the correlation between the tissue-specific expression and 
the target organism (antagonist) specificity of these molecules. In addition, conserved transcriptional circuitries 
involved in developing sexual organs in fungi are not characterized, despite the growing number of gene expression 
datasets available from reproductive and vegetative tissue. The aims of this study were: first, to evaluate the tissue 
specificity of defense gene expression in the model mushroom Coprinopsis cinerea and, second, to assess the degree of 
conservation in transcriptional regulation during sexual development in basidiomycetes. 

Results: In order to characterize the regulation in the expression of defense loci and the transcriptional circuitries 
controlling sexual reproduction in basidiomycetes, we sequenced the poly (A)-positive transcriptome of stage 1 
primordia and vegetative mycelium of C. cinerea A43mutB43mut. Our data show that many genes encoding predicted 
and already characterized defense proteins are differentially expressed in these tissues. The predicted specificity of 
these proteins with regard to target organisms suggests that their expression pattern correlates with the type of 
antagonists these tissues are confronted with. Accordingly, we show that the stage 1 primordium-specific protein 
CC1G_1 1805 is toxic to insects and nematodes. Comparison of our data to analogous data from Laccaria bicolor and 
Schizophyllum commune revealed that the transcriptional regulation of nearly 70 loci is conserved and probably 
subjected to stabilizing selection. A Velvet domain-containing protein was found to be up-regulated in all three fungi, 
providing preliminary evidence of a possible role of the Velvet protein family in sexual development of basidiomycetes. 
The PBS-soluble proteome of C cinerea primordia and mycelium was analyzed by shotgun LC-MS. This proteome data 
confirmed the presence of intracellular defense proteins in primordia. 

Conclusions: This study shows that the exposure of different tissues in fungi to different types of antagonists shapes 
the expression pattern of defense loci in a tissue-specific manner. Furthermore, we identify a transcriptional circuitry 
conserved among basidiomycetes during fruiting body formation that involves, amongst other transcription factors, 
the up-regulation of a Velvet domain-containing protein. 
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Background 

The last eukaryotic common ancestor (LECA) was facul- 
tatively sexual and evolved nearly 1.5 billion years ago in 
the Proterozoic eon [1]. Sexual reproduction (SR) shares 
common features across the eukaryotic lineage such as 
ploidy changes, meiotic recombination and cell-cell rec- 
ognition between gametes followed by cellular fusion 
and zygote formation [2]. An increase in the genetic di- 
versity of the population, making it more adaptable to 
changing environmental conditions, as well as the dilu- 
tion of deleterious mutations out of the gene pool are 
the most remarkable evolutionary innovations achieved 
by SR [3]. Despite all these clear benefits, sex is energet- 
ically expensive and entails a higher chance of genetic 
and organelle conflicts [2]. 

During sexual reproduction, fungi undergo dramatic 
morphological changes driven by environmental condi- 
tions such as light, nutrient availability and grazing by 
predators [4]. In basidiomycetes, mushroom develop- 
ment starts with intense localized hyphal branching 
leading to the formation of hyphal knots. These branch- 
ing hyphae further aggregate to form 1-2 mm secondary 
nodules where cell differentiation leads to the establish- 
ment of bipolar primordia containing all the tissues ob- 
served in the mature fruiting body [5]. As a last step, 
primordia develop to mature fruiting bodies mainly by 
cellular expansion [6]. Due to their hyphal density, prim- 
ordia and fruiting bodies are attractive to predators in- 
cluding mollusks, arthropods and nematodes [4]. 

Coprinopsis cinerea has been used as a model basidio- 
mycete since the mid-1950s [7] due to its saprobic life- 
style, its rapid growth and the feasibility of producing 
fruiting bodies under defined laboratory conditions [6]. 
In nature, C. cinerea grows on horse dung [6], a eu- 
trophic substrate rich in competing microorganisms, 
such as Firmicutes, Bacteroidetes and Proteobacteria 
[8]. The recent sequencing of the C. cinerea genome 
[9] allows the study of this organism gene expression 
on transcriptome and proteome level at different de- 
velopmental stages or under a variety of environmental 
settings. 

Morphological changes and environmental signals dur- 
ing fruiting body formation in C. cinerea are well de- 
scribed [6]; nonetheless, comparably little is known about 
the molecular machinery driving sexual reproduction pro- 
cesses in this basidiomycete. Recently, mutations blocking 
fruiting body development at different stages or altering 
mushroom morphology were identified [10-15]. In addi- 
tion, dst2 and dstl, encoding a blue-light photoreceptor 
and a flavin adenine dinucleotide-binding protein, were 
shown to play a role in blue light sensing. In agreement 
with previous experiments, strains defective in these two 
proteins were unable to form fruiting bodies, showing that 
blue light is an essential environmental trigger of 



mushroom development [16]. Most recently, C. cinerea 
strains carrying mutations in the putative component of 
the SWI/SNF chromatin remodeling complex sn/S 
(CC1G_15539) were shown to be defective in fruiting initi- 
ation, suggesting that epigenetic reprogramming of loci 
occurs during fruiting body formation [17]. 

Aerial fruiting bodies are an attractive prey for preda- 
tors and thus are protected by a battery of defense mole- 
cules (toxins) including proteins [18-21], peptides [22] 
and secondary metabolites [23-25]. Some of these toxins 
are known to be specifically produced in the fruiting 
body and not in the vegetative mycelium [26,27]. For in- 
stance, cytoplasmic lectins showing a broad range of 
non-self-carbohydrate specificities, also referred to as 
fruiting body lectins due to their specific expression pat- 
tern, have been shown to exert toxicity to nematodes, 
insect larvae and amoeba [18]. Protein-mediated inhib- 
ition of serine proteases [20], proteolytic degradation of 
predator-derived proteins [28,29] and sequestration of 
bio tin [19] are other strategies of basidiomycetous fruit- 
ing bodies to dissuade predators. Vegetative mycelium, 
in contrast, digests extracellular carbon macromolecules 
in the growth substrate, such as cellulose and lignin, into 
smaller degradation products which are absorbed by the 
growing hyphae. At the same time, these smaller mole- 
cules become available for competing bacteria that profit 
from the fungal enzymatic machinery [30]. As a re- 
sponse to bacterial competitors, fungi have evolved se- 
creted antimicrobial proteins. Mygind and colleagues 
presented recent evidence that vegetative mycelia of 
fungi secrete cysteine-stabilized antibacterial peptides 
which play a role in the arms race with competing 
bacteria [31]. 

Using RNA-seq in C. cinerea, we show evidence sug- 
gesting a role of the Velvet protein regulon in sexual de- 
velopment of basidiomycetes. Our data supports the 
existence of a conserved transcriptional circuitry in ba- 
sidiomycetes fruiting bodies consisting of at least 60 
orthologous genes probably involved in mushroom de- 
velopment and function. In addition, our data reveals 
the existence of two different sets of fungal defense pro- 
teins in vegetative and sexual organs matching the type 
of competitors and predators by which these structures 
are challenged in nature. The transcriptome data is sup- 
ported by the first partial shotgun mass spectrometry 
catalog of proteins present in C. cinerea stage 1 primor- 
dia (SIP) and vegetative mycelium (VM). 

Results 

Differential gene expression during fruiting body 
development in C. cinerea 

Four different cDNA libraries were sequenced with a 
final 7.32 Gb mapped data output. Approximately 95% 
of the open reading frames (ORFs) in the genome of 
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C. cinerea A43mutB43mut (AB) were transcribed using 
five reads/ORF as the minimal threshold for a locus to 
be considered as expressed (Table 1). Differential gene 
expression at the RNA level between SIP and VM was 
examined (Additional file 1: Table SI). Eleven percent of 
the annotated ORFs in the C. cinerea genome were 
found to be differentially expressed, 795 (6%) and 679 
loci (5%) in VM and SIP, respectively, using fold change 
8 and Fisher's exact test p-value < 0.05 as thresholds 
(Figure 1, Additional file 2: Table S2). The number of 
differentially transcribed loci in these two developmen- 
tal stages increased to 2522 in VM and 3209 in SIP 
when a fold change threshold of 2 was set, correspond- 
ing to approximately 45% annotated ORFs in the gen- 
ome of C. cinerea. 

To validate these data, qRT-PCR of four selected loci 
(CC1G_10318, CC1G_09480, CC1G_05299 and CC1G_ 
11805) was performed. Although S1P/VM expression ra- 
tios did not exactly match those observed by RNA-seq, 
qRT-PCR results showed the same trend of differential 
gene expression for these loci during C. cinerea develop- 
ment as observed by SOLiD RNA-seq (Additional file 3: 
Figure SI). In addition, we verified that the expression 
values of reported housekeeping loci commonly used in 
qRT-PCR normalization were in a range indicative of con- 
stitutive expression (-2 < log2(SlP/VM) < 2) (Additional 
file 4: Table S3). 

Enrichment of PSTBLAST-derived functional annota- 
tion terms was visualized using Wordle (© IBM Corpor- 
ation) [32] after excluding frequent non-informative 
terms. By far, the most commonly assigned annotation 
term was "hypothetical protein" in SIP (275 loci) and 
VM (426 loci); nonetheless, there was a significant 
enrichment of N-methyltransferase (histone-lysine N- 
methyltransferase) tags in SIP (7 loci). Enrichment of 
functional annotation terms also shows that different 
sets of cytochromes, kinases, dehydrogenases, trans- 
porters and hydrophobins are specifically expressed in 
vegetative mycelium or young fruiting bodies (Figure 1). 
Functional annotation clustering using DAVID Bioinfor- 
matics Resources 6.7 [33] confirmed the enrichment of 

Table 1 General features of the C. cinerea SI P and VM 



transcriptomes 



Library size (reads) 


Biological 


Biological 


Expression (%)* 


replicate 1 


replicate 2 




S1P 26243754 


25653679 


95.55 


VM 29374279 


23269631 


95.13 


Average SOLiD read size (b) 


70 




Sequenced material (Gb) 


7.32 





*Mean percentage of annotated ORFs detected to be expressed using a 
minimal threshold 5 reads/ORF. 



protein methyltransferases and hydrophobins in SIP with 
11 and 6 fold enrichment compared to the occurrence of 
these functional categories in the genome of C. cinerea 
(Additional file 5: Table S4). 

Tissue-specific expression of C. cinerea defense proteins 

Several cytoplasmic defense lectins and protease inhibi- 
tors such as C. cinerea lectin 2 (CCL2, CC1G_11781) 
[34], Coprinopsis galectin 1 (CGL1, CC1G_05003) [35] 
and two paralogous serine protease inhibitors from 
Coprinopsis (Cospin; CC1G_09479 and CC1G_09480) 
[20] were found to be highly up-regulated in C. cinerea 
SIP (Table 2). Out of these previously characterized 
genes coding for proteins with nematotoxic and insecti- 
cidal activity, Cospin and CCL2 were found to be among 
the top 50 most highly transcribed and differentially 
expressed loci in SIP. Among these, two loci, CC1G_ 
10318 and CC1G_11805, encoding homologous proteins 
with a predicted aerolysin/ETX pore-forming domain as 
well as loci encoding proteins with suspected antibacter- 
ial and antifungal function including two peptidoglycan 
binding proteins, the toxin component of a bacterial 
toxin-antitoxin system and a Thaumatin-like protein 
were found. In contrast to S IP-specific defense proteins, 
defense proteins specifically up-regulated in VM were 
mainly secreted. These proteins included three putative 
lysozymes, several proteins containing a CFEM domain 
(PF05730) whose structure resembles cysteine-stabilized 
antibacterial peptides, and two representatives of the 
cerato-platanin family of secreted proteins (Table 3). Lat- 
ter protein family is expanded in basidiomycetes and has 
recently been implicated in interactions of dikaryotic 
fungi with other organisms [36,37]. During sexual devel- 
opment on herbivore dung, C. cinerea is exposed to a 
succession of antagonists (predators and competitors) 
colonizing this substrate. The differential expression of 
cytoplasmic and secreted defense proteins in SIP and 
VM, respectively, might reflect the prevalent types of an- 
tagonists with which these tissues of C. cinerea are 
confronted. 

The S1 P-specific aerolysin/ETX pore-forming 
domain-containing protein CC1G_11805 is toxic to 
nematode and insect larvae 

CC1G_11805, CC1G_10318 and CC1G_08369 encode 
three homologous 30-40 kDa proteins containing a pre- 
dicted aerolysin/ETX pore-forming domain (Figure 2A). 
This domain is homologous to the one present in the in- 
secticidal epsilon toxin (ETX) from Clostridium perfrin- 
gens [38] and distantly related to aerolysin pore-forming 
toxins [39]. In order to test the significance of the up- 
regulation of CC1G_11805 and CC1G_10318 in SIP 
with regard to fungal defense, we cloned and recombi- 
nantly expressed CC1G_11805 in Escherichia coli and 
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Figure 1 Differential gene expression during early fruiting body development in C. cinerea. (A) Volcano plot illustrating the relative 
expression (log2(S1 P/VM)) and statistical significance of the 13342 loci in the genome of C. cinerea between VM and SI P. Fisher's exact test -loglO 
p-value > 1 .3 (p-value < 0.05) and log2(S1 P/VM) +/- 3 (8 fold change) were used as thresholds of differential gene expression. Grey boxes 
comprise genes significantly up-regulated in SIP (right) or VM (left). Tag clouds showing enriched PSI-BLAST functional annotation terms from 
up-regulated loci in S1P (B) and VM (C) were produced in Wordle (© IBM Corporation) after removing frequently appearing tags. 



Table 2 SI P-specific defense loci 



Locus 


Fold S1P/VM 


Functional annotation 


*p-value 


tSignalP 


tTMHMM 


CC1G_09480 


2426 


Cospinl 


1.26E-305 


N 


N 


CC1GJ1781 


1939 


CCL2 


4.64E-302 


N 


N 


CC1G_09479 


692 


Cospin2 


2.38E-132 


N 


N 


CC1GJ2219 


237 


Related to Velvet A protein 


6.42E-43 


N 


N 


CC1G_07937 


83 


Ricin B-fold protein 


1.08E-148 


Y 


N 


CC1G_1 1778 


69 


CCL1 


7.80E-137 


N 


N 


CC1GJ0318 


68 


Pore-forming protein 


8.08E-91 


N 


N 


CC1G_07956 


66 


Peptidoglycan-binding domain 1 protein 


1 .70E-48 


Y 


1 


CC1GJ4321 


35 


Hemolysin 


9.81 E-77 


N 


3 


CC1G_08484 


34 


Cercosporin toxin biosynthesis protein 


4.71E-115 


N 


N 


CC1G_06959 


31 


Thaumatin-like protein 


3.50E-79 


Y 


N 


CC1GJ3099 


15 


Peptidoglycan-binding domain 1 protein 


4.40E-17 


N 


N 


CC1G_05003 


12 


CGL1 


7.13E-37 


N 


N 


CC1G_1 1 123 


1 1 


Toxin-antitoxin system, toxin component 


1 .84E-33 


N 


N 


CC1GJ1246 


11 


Ricin B-fold protein 


7.70E-28 


Y 


N 


CC1GJ1805 


11 


Pore-forming protein 


1 .69E-40 


N 


N 



*Fisher's exact test. 

tN or Y indicates the lack or presence of a signal peptide. 

indicates a lack of transmembrane helices, while a number corresponds to the amount of transmembrane helices predicted in the ORF. 
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Table 3 VM-specific defense loci 



Locus 


Fold VM/S1P 


PSI Blast 


*p-value 


tSignalP 


+TMHMM 


CC1 G_05299 


2076 


Ricin B-fold protein 


2.95E-240 


N 


N 


CC1GJ0614 


1481 


CFEM domain-containing protein 


9.18E-166 


Y 


N 


CC1GJ5645 


851 


GFEM domain-containing protein 


1.51E-194 


Y 


N 


CC1G_09154 


639 


Cerato-platanin protein 


6.24E-233 


Y 


N 


CC1G_13813 


527 


CFEM domain-containing protein 


6.90E-195 


Y 


N 


CC1G_0S638 


67 


Peptidoglycan-binding domain 1 protein 


1.19E-119 


Y 


1 


CC1G_0S246 


42 


Ricin B-fold protein 


7.71E-125 


Y 


N 


CG1G_09155 


24 


Cerato-platanin protein 


3.27E-75 


Y 


N 


CG1G_09421 


21 


Terpenoid synthase 


5.98E-27 


N 


N 


CC1G_1 1847 


17 


Lysozyme 


2.00E-08 


Y 


N 


CC1G_08066 


16 


Ricin B-fold protein 


1.01 E-45 


Y 


N 


CC1G_08310 


13 


Lysozyme 


6.90E-09 


Y 


N 


CC1GJ5739 


10 


Ricin B-fold protein 


1.54E-59 


Y 


N 


CC1G_03046 


8 


Lysozyme 


3.68E-05 


Y 


N 



*Fisher's exact test. 

IN or Y indicates the lack or presence of a signal peptide. 

indicates a lack of transmembrane helices, while a number corresponds to the amount of transmembrane helices predicted in the ORF. 



assessed the toxicity of the protein by feeding the re- 
combinant bacteria to nematode and insect larvae as de- 
scribed previously [40]. Results in Figure 2B show that 
CC1G_11805 was expressed in soluble form in E. coli. 
Feeding of CClG_11805-expressing bacteria to L2 larvae 
of the mosquito Aedes aegypti lead to their death after 
96 h (Figure 2C). Vector control-containing and fungal 
lectin-expressing E. coli were used as negative and posi- 
tive controls, respectively, in these experiment [18,29]. 
Similarly, feeding of CClG_11805-expressing bacteria 
was found to significantly impair larval development of 
the C. elegans N2 wildtype and pmk-l(km25) mutant 
strains (Figure 2D); in agreement with the previously re- 
ported higher susceptibility of latter strain to different 
kinds of abiotic and biotic stresses including other 
nematotoxic fungal defense proteins [35], CC1G_11805 
was shown to be more toxic to C. elegans pmk-l(km25) 
than to N2. These results support previous observations 
[18] showing that the expression of defense proteins di- 
rected against nematodes and insects is significantly in- 
creased in C. cinerea sexual organs. 



Aspergillus nidulans [41,42]; therefore, we took a closer 
look at the expression pattern of C. cinerea orthologs of 
the Velvet protein family and the proteins regulated by it 
[41,42]. CC1G_12219, a Velvet domain-containing pro- 
tein bearing an NLS is highly induced in SIP (S1P/VM = 
237). In total, three out of six Velvet domain-containing 
proteins encoded in the C. cinerea genome (CC1G_12219, 
CC1G_06962 and CC1G_14883) were found to be devel- 
opmentally regulated (Additional file 7: Figure S3). Among 
the other members of the Velvet protein regulon, two ho- 
mologs of the VeA-regulated activator NsdD (CC1G_ 
06391 and CC1G_06265) were also up-regulated in SIP 
(S1P/VM = 14 and 2.4, respectively). Moreover, the C 
cinerea homeodomain transcription factor STE-12 (CC1G_ 
02207) showed a moderate up-regulation in primordia 
(S1P/VM = 2). Finally, CC1G_07060, an ortholog of the A. 
nidulans repressor of sexual development rosA, was 34 
fold up-regulated in VM compared to SIP. Taken to- 
gether, the differential expression of several members of 
the Velvet protein regulon suggests that these genes may 
play a role in sexual development of basidiomycetes. 



Differential expression of loci involved in sexual 
development of filamentous fungi 

Various genes were recently described as playing a role 
at different stages of fruiting body formation of C. 
cinerea [10-17]. With the exception of the cyclopropane 
fatty acid synthase cfsl, none of these loci were found to 
be differentially transcribed during C. cinerea sexual de- 
velopment (Table 4, Additional file 6: Figure S2). The 
Velvet protein regulon plays a major role in the control 
of sexual vs. asexual development in the ascomycete 



Conserved transcriptional circuitry during fruiting body 
development among basidiomycetous fungi 

In order to explore a hypothetical conserved gene ex- 
pression circuitry during the formation of SR structures 
in basidiomycetes, a comparative transcriptomic analysis 
by hierarchical clustering including C. cinerea, S. com- 
mune [43] and L. bicolor [44] was performed. Ortholo- 
gous genes in S. commune and L. bicolor were assigned 
from a PSI-BLAST search (hits with the lowest E- 
value < 0.005) to 446 and 378 up- and down-regulated 
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Figure 2 The S1P-specific aerolysin/ETX pore-forming domain-containing protein CC1G_1 1805 is toxic for C. elegans and A. aegypti. 

(A) Multiple sequence alignment of the aerolysin/ETX pore-forming domain-containing proteins C. perfringens epsilon toxin (Cper_etx), CC1G_1 1805, 
CC1GJ0308 and CC1G_08369. (B) Expression of soluble CC1GJ 1805 in soluble form in £ coli BL21. Whole cell protein extracts in PBS were produced 
from non-induced (WCE) and induced (WCE + IPTG) £ coli BL21/pET24-CC1G_1 1805 cultures. Supernatants from consecutive low (LS) and high (HS) 
speed centrifugations were collected and run on a 12% SDS-PAGE to evaluate CC1G_1 1805 expression and solubility. CC1G_1 1805-mediated 
toxicity for A. aegypti (C) and C. elegans (D) larvae was assessed as described previously [40]. £ coli BL21 expressing the previously characterized fungal 
lectin CGL2 was used as positive control for toxicity against A. aegypti Rockefeller. IPTG-induced £ coli bearing an empty pET24 vector (EV) was used 
as a negative control. Columns represent the mean percentage of surviving insect larvae (C), or worms either dying or reaching the indicated 
developmental stage (D) from 4 replicates each. SDs are shown as error bars. Dunn's multiple comparisons were used to test the statistical significance 
of the toxicity observed in the A. aegypti assay. Mann Whitney test was used to compare the percentage of worms reaching each developmental stage 
when treated with EV or CC1GJ 1805. *: 0.01 < p-value < 0.05; **: 0.001 < p-value < 0.01. 



loci in C. cinerea SIP, respectively. Massively parallel 
signature sequencing tags/million from S. commune 
monokaryotic mycelium and stage 1 primordia (ScSIP); 
as well as L. bicolor free-living monokaryotic mycelium 
(2 replicates) and young fruiting bodies (L&YFB, one rep- 
licate) robust multichip average values were retrieved 



from original studies [43,44]. log2 expression ratios were 
calculated using monokaryotic mycelium expression 
values from these conserved loci as denominators and a 
centroid linkage hierarchical clustering analysis was 
computed (Figure 3; Additional file 8: Table S5). Four 
clusters corresponding to 37 up-regulated loci in early 
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Table 4 Most genes necessary for fruiting body formation in C. cinerea are not up-regulated in SIP 



Locus 


Fold S1P/VM 


Gene 


*p-value 


Mutant phenotype 


Reference 


CC1 GJB287 


0.8 


rmtl 


1.43E-01 


No clamp connections 


Nakazawa [15] 


CC1G_00975 


1.4 


ubc2 


1.75E-02 


No hyphal knot formation 


Nakazawa 201 1 


CC1G 15539 


1.1 


snf5 


3.98E-01 


B ocked at hyphal knot 


Ando [1 7] 


CC1GJ1387 


64.0 


cfsl 


5.25E-142 


Blocked at initials 


Liu [10] 


CC1G_06825 


1.2 


dst2 


9.66E-02 


Blocked at stage 2 primordia 


Kamada [16] 


CC1G_08609 


2.6 


dstl 


1.21E-13 


Blocked at stage 2 primordia 


Kamada [16] 


CC1GJ0193 


1.5 


eln3 


1 .40E-03 


Blocked at immature mushroom 


Arima [11] 


CC1G_04713 


2.3 


eln3 


5.18E-08 


Blocked at immature mushroom 


Arima [11] 


CC1G_06451 


2.5 


eln3 


9.03 E- 13 


Blocked at immature mushroom 


Arima [11] 


CG1G_01334 


0.8 


expl 


3.65E-02 


Blocked before mushroom decay 


Muraguchi [12] 



•Fisher's exact test. 



sexual development in basidiomycetes were identified 
(Pearson correlation coefficients: 0.95, 0.99, 0.99 and 
0.99, respectively). In addition, two different gene groups 
(29 loci) were found to be down-regulated in early stages 
of fruiting body formation when compared to vegetative 



mycelium of the three species (Pearson correlation coef- 
ficients: 0.99 and 0.94, respectively) (Table 5). 

This hierarchical clustering approach of gene expres- 
sion in orthologous loci during sexual development 
showed that loci involved in core meiotic functions, such 



ScS1P CcS1P LbYFB 



B 



ScS1P CcS1P ZJbYFB 








log2(fold) log2(fold) 



-5> 0 5< -5> 0 5< 

Figure 3 Conserved expression pattern of some gene clusters during fruiting body formation among three different basidiomycetous 
species. Comparative transcriptome analysis of early sexual development in 5. commune 4-39/4-40, C. cinerea AB and L bicolor S238N-H82 was performed. 
Orthologous genes (best PSI-BLAST hits with E-values < 0.005) to only 446 and 378 up- and down-regulated loci in SI P, respectively, were found in 
5. commune 4-39/4-40 and L bicolor S238N-H82. (A) Fruiting body transcriptome data were retrieved in the form of log2 ratios for these loci from the raw 
data sets and analyzed by centroid-linkage hierarchical clustering. The red sidelines mark the presence of 6 loci clusters which are consistently up (1-4)- or 
down (5 and 6)-regulated in young fruiting bodies (YFB) and stage 1 primordia (S1P) of the three species compared. (B) Close-up view of the up- and 
down-regulated clusters highlighted by red bars in (A). Orthologous genes in C. cinerea, 5. commune and L bicolor corresponding to the numbered 
clusters, as well as the functional annotation of these loci and the cluster-associated Pearson correlation coefficients can be found in Table 5. 
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Table 5 Up- and down-regulation of orthologous loci is conserved in stage 1 primordia and young fruiting bodies of 
three basidiomycete species 



Cluster 


C. cinerea locus 


S. commune locus 


L. bicolor locus 


Functional annotation* 


1(PCt: 0.95) 


CC1G_01874 


SCHCODRAFT_73392 


LACBIDRAFT_301146 


Hypothetical protein 




CC1G_02466 


SCHCODRAFT_1 09665 


LACBIDRAFT_297494 


Hypothetical protein 




CC1G_00973 


SCHCODRAFT_83399 


LACBIDRAFT_293318 


Chitin deacetylase 




CC1GJ0475 


SCHCODRAFT_57566 


LACBIDRAFT_295571 


Aromatic peroxygenase precursor 




CC1GJ4095 


SCHCODRAFT_48101 


LACBIDRAFT_294461 


Hypothetical protein 




CC1G_03094 


SCHCODRAFT_52606 


LACBIDRAFT_315775 


Hypothetical protein 




CC1G_06074 


SCHCODRAFT_76530 


LACBIDRAFT_293318 


Carbohydrate esterase family 4 protein 




CC1GJ0044 


SCHCODRAFT_1 08884 


LACBIDRAFT_306386 


Related to S.pombe pac2 protein 




CC1G_00753 


SCHCODRAFT_234371 


LACBIDRAFT_305130 


Hypothetical protein 




CC1GJ4786 


SCHCODRAFT_52306 


LACBIDRAFT_246776 


MFS nicotinic acid transporter Tna1 




CC1G_01301 


SCHCODRAFT_72462 


LACBIDRAFT_309640 


Hypothetical protein 




CC1G_05059 


SCHCODRAFT_52448 


LACBIDRAFT_308735 


Symbiosis-related protein 




CC1G_1 1437 


SCHCODRAFT_67374 


LACBIDRAFT_317173 


Putative aquaporin 6 




CC1G_08178 


SCHCODRAFT_237423 


LACBIDRAFT_323328 


Thioesterase family protein 




CC1G_08180 


SCHCODRAFT_237423 


LACBIDRAFT_323328 


Thioesterase 




CC1G_01879 


SCHCODRAFT_257636 


LACBIDRAFT_326397 


mei2 protein Piriformospora 




CC1G_05223 


SCHCODRAFT_70203 


LACBIDRAFT_250946 


DUF1275 domain protein 




CC1G_12219 


SCHCODRAFT_28806 


LACBIDRAFT_317102 


Related to velvet A protein 




CC1G_10750 


SCHCODRAFT_86141 


LACBIDRAFT_1 84665 


Glutathione S-transferase 




CC1G_04060 


SCHCODRAFT_13677 


LACBIDRAFT_1 80892 


Hydrophobin-251 




CC1G_10471 


SCHCODRAFT_232646 


LACBIDRAFT_295571 


Aromatic peroxygenase precursor 




CC1G_.12330 


SCHCODRAFT_16019 


LACBIDRAFT_256021 


S-layer domain-containing protein 


2 (PC: 0.99) 


CC1G_00069 


SCHCODRAFT_48115 


LACBIDRAFT_312065 


C factor cell-cell signaling protein 




CC1G_00217 


SCHCODRAFT_46102 


LACBIDRAFT_323571 


Non-Catalytic module family EXPN protein 




CC1G_03515 


SCHCODRAFT_82058 


LACBIDRAFT_328112 


Hydrophobin-like protein 




CC1G_04659 


SCHCODRAFT_1 10821 


LACBIDRAFT_298314 


Hypothetical protein 




CC1G_05652 


SCHCODRAFT_232455 


LACBIDRAFT_309898 


TKL/TKL-ccin protein kinase 




CC1G_09209 


SCHCODRAFTJ 10821 


LACBIDRAFT_298284 


BTB/POZ domain containing protein 




CC1G_10966 


SCHCODRAFT_258323 


LACBIDRAFT_236299 


rCop c3 


3 (PC: 0.99) 


CC1G_02571 


SCHCODRAFT_76530 


LACBIDRAFT_293318 


Chitin deacetylase 




CC1G_03037 


SCHCODRAFT_53388 


LACBIDRAFT_239714 


Monocarboxylate transporter 




CC1G_03320 


SCHCODRAFT_39718 


LACBIDRAFTJ45110 


Spoil 




CC1G_09616 


SCHCODRAFT_235946 


LACBIDRAFT_310809 


Prenyl cysteine carboxyl methyltransferase 




CC1G_06563 


SCHCODRAFT_46720 


LACBIDRAFT_243581 


Exo-beta-1,3-glucanase 




CC1GJ2688 


SCHCODRAFT_55636 


LACBIDRAFT_296585 


Aldo/keto reductase 


4 (PC: 0.99) 


CC1G_03949 


SCHCODRAFT_258034 


LACBIDRAFT_298461 


Salicylate hydroxylase 




CC1G_05817 


SCHCODRAFT_258034 


LACBIDRAFT_298461 


Salicylate 1-monooxygenase 


5 (PC: 0.99) 


CC1G_06868 


SCHCODRAFT_53362 


LACBIDRAFT_1 82606 


Endo-1 ,3(4)-beta-glucanase 




CC1G_07550 


SCHCODRAFT_56140 


LACBIDRAFT_291657 


Lipase 




CC1G_05781 


SCHCODRAFT_72461 


LACBIDRAFT_291413 


OrdA protein 




CC1G_08860 


SCHCODR AFT_80 1 8 1 


LACBIDRAFT_1 72283 


Cation/H + exchanger 




CC1G_09431 


SCHCODRAFT_49922 


LACBIDRAFTJ 90903 


O-methylsterigmatocystin oxidoreductase 




CC1G_11000 


SCHCODRAFT_1 07054 


LACBIDRAFT_310382 


Hypothetical protein 
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Table 5 Up- and down-regulation of orthologous loci is conserved in stage 1 primordia and young fruiting bodies of 
three basidiomycete species (Continued) 



6 (PC: 0.94) 



cc 


IG. 


02440 


SCHCODRAFT 


.103563 


LACBIDRAFT 


.308947 


ras GEF 


cc 


IG. 


07544 


SCHCODRAFT 


.80832 


LACBIDRAFT 


.243185 


DUF89-containing protein 


cc 


IG. 


05410 


SCHCODRAFT. 


.104205 


LACBIDRAFT 


.301153 


Hypothetical protein 


cc 


IG. 


.02628 


SCHCODRAFT. 


.82883 


LACBIDRAFT 


.327335 


C2h2 conidiation transcription factor FlbC 


cc 


IG. 


.1 1 894 


SCHCODRAFT. 


.68168 


LACBIDRAFT 


.296675 


basl , putative 


cc 


IG. 


05329 


SCHCODRAFT. 


.85265 


LACBIDRAFT 


.300118 


Hypothetical protein 


cc 


IG. 


.09061 


SCHCODRAFT. 


.27314 


LACBIDRAFT 


.308057 


MNNG and nitrosoguanidine resistance protein 


cc 


IG. 


.02380 


SCHCODRAFT. 


.77933 


LACBIDRAFT 


.314347 


Clavaminate synthase-like protein 


cc 


IG. 


05793 


SCHCODRAFT. 


.85265 


LACBIDRAFT 


.300118 


Hypothetical protein 


cc 


IG. 


.14014 


SCHCODRAFT. 


.49922 


LACBIDRAFT 


.248655 


O-methylsterigmatocystin oxidoreductase 


cc 


G_ 


.11136 


SCHCODRAFT. 


.232299 


LACBIDRAFT 


.294457 


Hypothetical protein 


cc 


IG. 


.00780 


SCHCODRAFT. 


.1 06836 


LACBIDRAFT 


.305061 


BGP, partial 


cc 


IG. 


03120 


SCHCODRAFT. 


.83385 


LACBIDRAFT 


.307143 


Endoglucanase 1 


cc 


IG. 


02257 


SCHCODRAFT. 


.46132 


LACBIDRAFT 


.295167 


Metalloprotease 


cc 


IG. 


04116 


SCHCODRAFT. 


.75425 


LACBIDRAFT 


.1 79508 


Glutathione S-transferase Gst3 


cc 


IG. 


.11081 


SCHCODRAFT. 


.75425 


LACBIDRAFT 


.1 79508 


Glutathione S-transferase 


cc 


IG. 


.12752 


SCHCODRAFT. 


.50449 


LACBIDRAFT 


.253329 


Cytochrome P450 


cc 


IG. 


J 2253 


SCHCODRAFT. 


.269973 


LACBIDRAFT 


.296037 


Nrgl -like Zn-finger transcription factor 


cc 


IG. 


.1 1 786 


SCHCODRAFT. 


.58473 


LACBIDRAFT 


.315853 


Hypothetical protein 


cc 


IG. 


.09938 


SCHCODRAFT. 


.80558 


LACBIDRAFT 


.330070 


Ferritin/ribonucleotide reductase-like-protein 


cc 


IG. 


.02724 


SCHCODRAFT. 


.255861 


LACBIDRAFT 


.295183 


Hypothetical protein 


cc 


IG. 


.05802 


SCHCODRAFT. 


.65657 


LACBIDRAFT 


.256254 


lndole-3-acetaldehyde dehydrogenase 


cc 


IG. 


.15739 


SCHCODRAFT. 


.76887 


LACBIDRAFT 


.146952 


Ricin B fold protein 



*Functional annotation as determined by PSI-BLAST. 
tPearson correlation coefficient. 



as mei2 (CC1G_01879) and spoil (CC1G_03320), are 
up-regulated during the early stages of fruiting body for- 
mation in basidiomycetes. 

Given the up-regulation of a conserved Velvet-domain 
containing protein in primordia of the analyzed basidio- 
mycetes, the conservation and expression pattern of 
Velvet-interacting proteins previously characterized in 
Aspergillus [42] was examined. With the exception of 
cry A, orthologs to these Velvet- interacting proteins can be 
found in C. cinerea, L. bicolor and S. commune (Additional 
file 9). Contrary to the conserved transcriptional reg- 
ulation observed for genes encoding Velvet domain- 
containing proteins, little conservation in expression is 
evident for the orthologs to rosAlnosA, stuA, nsdD, ppoA, 
laeA or jphA (Additional file 10: Figure S4). Nevertheless, 
conserved down-regulation during sexual development is 
shown for velB and kapA, suggesting that a fraction of the 
Velvet-associated regulon described for ascomycetes 
might be playing a role during sexual development in ba- 
sidiomycetes (Additional file 10: Figure S4). Similarly, 
down-regulation of transcription factors such as flbC (a 
paralog to C2H2 known to regulate development in A. 



nidulans [45]), nrgl and basl during early mushroom for- 
mation was conserved between the three basidiomycetous 
species. 

Transcription factor mutants (Bril, Homl, Gatl, Fst3, 
C2h2, Fst4 and Hom2) altering normal sexual develop- 
ment in S. commune have been described [46,47]. Al- 
though orthologs of these transcription factors are 
present in C. cinerea and L. bicolor, the respective genes 
do not show a conserved expression pattern, indicating 
that there might be differences in the interplay between 
these factors during fruiting body formation between 
different basidiomycetes (Additional file 11: Figure S5). 

Shotgun MS analysis of PBS-soluble VM and SIP proteins 

In order to confirm protein expression of some of the 
genes found to be differentially expressed by RNA-seq, 
we assessed the PBS-soluble C. cinerea proteome in VM 
and SIP by LC-MS. The analysis detected peptides cor- 
responding to a total of 493 proteins in the samples, in- 
cluding 41, 141 and 311 proteins in SIP, VM or both, 
respectively (Additional file 12: Table S6). This analysis 
is highly biased towards abundant and soluble proteins 
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and is likely to have failed to detect peptides from most 
of the loci identified by RNA-seq. Nevertheless, this 
method allowed us to confirm the translation of 12 tran- 
scripts up-regulated in SIP and 50 transcripts up- 
regulated in VM into proteins in the respective tissues 
(Additional file 13: Table S7). 

Discussion 

Our data shows a transcriptional switch during the dif- 
ferentiation of primordia in C. cinerea comprising 11% 
of the protein-encoding genome being up- or down- 
regulated. Differential transcription analyses carried out 
in other basidiomycetous fungi, such as Agrocybe aeger- 
ita, Cordyceps militaris and Ganoderma lucidum, have 
revealed similar or even more extensive transcriptional 
switches [48-50]. Taking the differential expression 
threshold used in the present study (log2 (fruiting body/ 
vegetative mycelium) > 3), 25% and 30% out of 18474 A. 
aegerita loci are up-regulated in fruiting bodies and 
vegetative mycelium, respectively [48]. Similarly, C. mili- 
taris up-regulates 40% loci during fruiting body forma- 
tion [49], while G. lucidum boosts the transcription of at 
least 27% measured loci in its sexual organs [50]. 

The enrichment among the up-regulated loci in SIP of 
the functional annotation tag "N-methyltransferase" (his- 
tone-lysine N-methyltransferase, 7 loci in total) suggests 
that epigenetic regulation of gene expression or chromo- 
some dynamics may play an important role during 
sexual development in C. cinerea. In Saccharomyces cer- 
evisiae, histone H3 lysine 4 trimethylation marks the 
sites where double-strand breaks, the first events during 
interhomolog recombination, occur [51]. Up-regulation 
of histone-lysine N-methyltransferases in SIP is in ac- 
cordance with karyogamy and early meiotic phases 
taking place at stage 1 primordia during fruiting body 
formation [6]. 

Several previously characterized nematotoxic and in- 
secticidal proteins [18,20] were found to be up-regulated 
in SIP. Among these proteins, lectins appear to be most 
abundant. The fruiting body-specific ricin B-like lectins 
CCL1 and CCL2 were previously shown to be toxic to 
C. elegans due to the binding to al,3-fucosylated N- 
glycan cores in the intestine of LI worms [34]. In ac- 
cordance to a previous study [52], we found that the 
tetrameric galectin CGL1 (CC1G_05003) was highly in- 
duced in C. cinerea fruiting bodies. Similar to CCL1 and 
CCL2, CGL1 and its isogalectin CGL2 (CC1G_05005) 
showed toxicity against C. elegans which was dependent 
on binding to a Gal[31,4Fucal,6-epitope on N-glycan 
cores present in the worm intestine [35]. Protease inhib- 
ition is another defense strategy against predation in C. 
cinerea fruiting bodies. We found the locus encoding the 
serine protease inhibitor Cospinl (CC1G_09480) and its 
isoprotein Cospin2 (CC1G_09479) among the top 50 



most highly transcribed and differentially expressed loci 
in SIP. Cospinl was previously shown to be toxic 
against Drosophila melanogaster larvae indicative of its 
role in fruiting body defense against arthropod predation 
[20]. In this study, we demonstrate that CC1G_11805 is 
toxic for C. elegans and A. aegypti larvae, analogous to 
the toxicity of aerolysin/ETX pore-forming domain- 
containing Bacillus sphaericus protein towards larvae of 
the mosquito Culex quinquefasciatus [53], suggesting 
that pore formation is another defense strategy of C. 
cinerea primordia against predators. The high absolute 
transcription of these defense loci in SIP can be ex- 
plained by the extensive resource allocation to these 
organs and the significance of these organs for repro- 
duction of the fungus. The expression of putative se- 
creted and antibacterial proteins in C. cinerea vegetative 
mycelium, on the other hand, probably reflects the con- 
frontation of this tissue with bacterial competitors al- 
though the antibacterial activity of these proteins has 
still to be demonstrated. In summary, the differential ex- 
pression of defense proteins in the different tissues of C. 
cinerea is an adaptation of the fungus to the different 
environmental challenges with which these tissues are 
confronted. 

Surprisingly, with the exception of cfsl, none of the 
genes previously shown to play a role in C. cinerea fruit- 
ing body formation were differentially expressed in stage 
1 primordia. Mutations in rmtl, ubc2 and snfS [14,15,17] 
block development before initials are formed (48 h before 
SIP develop). In contrast, dstl, dst2, eln3 and expl are in- 
volved in processes taking place in stage 2 primordia (24 h 
after SIP), immature fruiting bodies (48-72 h after SIP) 
and decaying fruiting bodies (72 h after SIP) [6,11,12,16]. 
Thus, a possible explanation of our results is that these 
genes might be regulated before or after the formation of 
SIP takes place. Intriguingly, the mutation of cfsl blocks 
development right at the transition between initials and 
SIP. This observed up-regulation of cfsl in SIP supports a 
function for this cyclopropane fatty acid synthase in the 
development of SIP from initials. 

The Velvet protein regulon, including genes nsdD, 
rosA, veA and stuA among others, coordinates sexual de- 
velopment and secondary metabolism in filamentous as- 
comycetes [41,42]. Velvet domains structurally resemble 
the RHD-like fold present in the transcription factor 
NF-k(3 that plays a central role in animal immunity, sug- 
gesting a common evolutionary origin for these two pro- 
tein families [54]. Overexpression of veA and nsdD in A. 
nidulans induces the formation of nursing Hiille cells 
surrounding the cleistothecia [41]. Homologous genes in 
C. cinerea, CC1G_12219 and CC1G_06391, showed high 
expression in SIP suggesting a function of these loci in 
the gene circuitry involved in fruiting body development 
in this fungus. On the other hand, rosA, an A. nidulans 
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transcription factor inhibiting sexual development in 
low-carbon culture, is expressed in A. nidulans asexual 
hyphae where it represses the transcription of sexual de- 
velopment regulators such as nsdD, veA and stuA [55]. 
Similarly, rasA-homologous C. cinerea gene CC1G_ 
07059, is expressed in vegetative mycelium where it 
might be inhibiting sexual development. Taken together, 
these results suggest a conserved role of the Velvet pro- 
tein regulon during sexual development in the asco- 
mycete A. nidulans and the basidiomycete C. cinerea. 

Genes encoding Velvet domain-containing proteins 
were also found to be up-regulated in fruiting bodies 
of S. commune and L. bicolor [43,44]. However, (with 
the exception of velB and kapA) genes encoding Velvet- 
associated proteins such as RosA/NosA, StuA, NsdD, 
PpoA, LaeA, FphA or CryA [42], do not show inter- 
species conservation of differential expression. Similar 
lack of transcriptional conservation between the three 
analyzed basidiomycete species was observed for ortho- 
logs to transcription factors shown to be involved in S. 
commune sexual development [46,47]. Taken together, 
this evidence suggests that basidiomycetes show species- 
specific divergence of transcriptional regulation in ortho- 
logous genes similar to plants [56]. 

A broader comparison of the transcriptomes of the 
three analyzed basidiomycete species, comprising all the 
genes differentially expressed in C. cinerea during sexual 
development, revealed the presence of a conserved gene 
regulation circuitry among basidiomycetes during fruit- 
ing body formation. As previously observed in fruiting 
bodies of ascomycetes [57], our orthology analysis 
showed that a large fraction of genes differentially 
expressed in SIP corresponds to loci which are not 
present in the basidiomycetes L. bicolor and S. commune. 
These results are in agreement with previous observa- 
tions in plants and animals that genes associated with 
sexual reproduction rapidly evolve [58,59]. The existence 
of clusters comprising conserved up- or down-regulated 
loci with little inter-species expression variability, sug- 
gests that regulation of these loci evolved under stabiliz- 
ing selection [60]. Sequence and expression conservation 
might imply an essential role of these genes in fruiting 
body development and sexual reproduction in basidio- 
mycetes. Interestingly, L. bicolor and C. cinerea ortho- 
logs to transcription factors described previously as 
important for sexual development in 5. commune, such 
as Bril, Homl, Gatl, Fst3, C2h2, Fst4 or Hom2 [46,47], 
do not show conserved transcriptional conservation in 
basidiomycetes, suggesting that novel regulatory path- 
ways related to sexual development did evolve once spe- 
ciation occurred. As a proof of concept, Traeger and 
collaborators found that the fruiting body specific tran- 
scription factor pro44 (orthologous to the velvet regulon- 
associated transcription factor nsdD) in the ascomycetes 



Sordaria macrospora and Pyronema confluens is a core 
regulator of perithecia maturation. S. macrospora deficient 
in pro44 was shown to be sterile and unable to produce 
mature perithecia [57]. In addition to a Velvet domain- 
containing protein, the expression of the RNA binding 
protein mei2 (CC1G_01879), a master meiosis regulator 
in yeast and plants [61,62] protecting meiosis-specific 
transcripts from degradation by the DSR-Mmi system, 
was increased in SIP or YFB of C. cinerea, S. commune 
and L. bicolor. Similarly, the transcriptional regulation of 
spoil, encoding a protein inducing meiotic recombination 
in S. cerevisiae and C. cinerea [63], was also conserved. 
This induction of meiosis regulators in multiple species 
reflects the role of mushrooms in the production and dis- 
persal of basidiospores. 

With regard to the comparative transcriptome analysis 
of the three basidiomycete species, it should be noted 
that in case of S. commune and L. bicolor, only data of 
monokaryotic vegetative mycelia was available [43,44], 
whereas in case of C. cinerea, the transcriptome of an 
isogenic homodikaryotic mycelium was determined. 
Thus, the degree of conserved transcriptional regulation 
during sexual development between C. cinerea, L. bicolor 
and 5. commune is potentially larger than observed. 

Lastly, the proteome of an organism provides a more 
direct image of its phenotype than the transcriptome 
[64]. Detection of proteins using label-free shotgun mass 
spectrometry fails to detect low abundance proteins in 
complex total extracts and allows only semi-quantitative 
estimation of relative protein amounts. These properties 
are in contrast to the superior standardization and sensi- 
tivity achieved by state of the art nucleotide sequencing 
technologies. Nevertheless, LC-MS spectra showed the 
presence of peptides derived from the nematotoxic lectin 
CCL1 (CC1G_11778) and the nematotoxic/insecticidal 
aerolysin/ETX pore-forming domain-containing protein 
CC1G_11805 in the PBS-soluble protein extract of C. 
cinerea stage 1 primordia, indicating that these cytoplas- 
mic toxins are indeed expressed at protein level in these 
organs. 

The recent development of gene targeting tools in C 
cinerea [13,65], will help to test whether some of the 
genes identified as differentially regulated during sexual 
development in multiple basidiomycete species, play a 
role in this process. 

Conclusions 

In this work, we show that sexual reproduction in C. 
cinerea A43mutB43mut involves the differential tran- 
scription of at least 11% of its protein-coding genome. 
Differentially transcribed genes include several genes 
coding for defense proteins that protect fruiting struc- 
tures and vegetative mycelia from predators and compet- 
itors, respectively. Moreover, our data infers a role of the 
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Velvet protein family during fruiting body formation 
in basidiomycetes and thus, a conserved role of this pro- 
tein family during sexual development in dikaryotic fila- 
mentous fungi. Finally, the result of the comparative 
transcriptome analysis of C. cinerea, S. commune and L. 
bicolor suggest that a conserved set of orthologous genes 
regulates sexual development in the phylum Basidiomy- 
cota. Additional experiments addressing the function of 
the differentially expressed gene products are required 
to confirm these hypotheses. 

Methods 

Strains and culture conditions 

The dikaryotic, self-compatible C. cinerea strain A43mut 
B43mut (AB) [66] was grown on 30 mL YMG plates 
(0.4% yeast extract, 1% malt extract, 50 mM glucose and 
1.5% agar) at 37°C in the dark for 96 h and transferred 
to 25°C, 90% humidity and 12 h photoperiod for fruiting 
body production. 1 to 2 mm AB SIP [6] were harvested 
after 72 h, flash-frozen in liquid nitrogen and stored at 
-80°C for later use. C. cinerea AB VM was grown in du- 
plicate on YMG plates covered with cellophane discs for 
96 h at 37°C in the dark and harvested independently 
before being flash frozen and stored at -80°C. 

Total RNA extraction 

Primordia and mycelia were lyophilized and SIPs were 
separated in two pools of 20 mg each. From each VM 
replicate or SIP pool, 20 mg dry material were lysed in 
three FastPrep FP120 homogenization steps of 45 s at 
4.5, 5.5 and 6.5 m/s in the presence of 250 mg 0.5 mm 
glass beads, cooling the samples for 5 min on ice be- 
tween steps. RNA was extracted using 1 mL Qiazol 
(Qiagen) and 0.2 mL chloroform ReagentPlus (Sigma- 
Aldrich). The solution was centrifuged at 12000 x g for 
15 min at 4°C; thereafter, RNA from the resultant aque- 
ous phase was washed on-column using the RNeasy 
Lipid Tissue Mini Kit (Qiagen) and eluted in 60 uL 
RNase-free water. Concentration and integrity of the 
purified RNA were determined with a Qubit (1.0) 
fluorometer (Life Technologies) and a Bioanalyzer 2100 
(Agilent), respectively. Samples with a 260/280 nm ratio 
of 1.8-2.1 and a 28S/18S ratio of 1.5-2 were later used 
in library construction. 

SOLiD 4 library construction 

Whole transcriptome libraries from two SIP pools and 
two VM replicates were produced using MicroPolyA 
Purist Kit (Ambion) and SOLiD Total RNA-Seq kit (Ap- 
plied Biosystems). Briefly, approximately 200 ng/sample 
poly (A)-positive RNA was enriched using MicroPoly 
(A) Purist Kit from 15-20 ug total RNA. Quality and 
concentration of the extracted poly (A)-positive RNA 
was re-assessed as described above, and poly (A) -positive 



RNA was digested with RNase III. Ligation of the 
adaptor mix and reverse transcription was performed 
following the manufacturer instructions. cDNA libraries 
were size selected for 150-250 bp fragments, amplified 
in 15-18 PCR cycles using barcoded adaptor primers 
and purified with PureLink PCR micro kit (Invitrogen). 
These barcoded cDNA libraries were then amplified by 
emulsion PCR from 0.5 pM template. Sequencing beads 
from the barcoded libraries were pooled and loaded on a 
SOLiD 4 slide (Applied Biosystems), according to manu- 
facturer's instructions. SOLiD ToP Sequencing chemistry 
was used to produce pair end (50 bp + 35 bp) sequencing 
reads. 

qRT-PCR validation 

RNA-seq results were validated by qRT-PCR. Single- 
stranded cDNA from one biological replicate per sample 
was synthesized using Transcriptor Universal cDNA 
Master (Roche) from 2 ug total RNA. 20 uL qRT-PCR 
reactions were mixed in three technical replicates per 
primer set and sample, containing 900 nM forward and 
reverse primers designed to span exon-exon junctions 
(Additional file 14: Table S8), 10 uL 2x FastStart Univer- 
sal SYBR Green Master (Rox, Roche) and 1 ng/uL cDNA 
template. qRT-PCR was performed in a Rotor-Gene 
3000 (Corbett Life Science) with the following thermal 
profile: a hold step at 95°C for 15 min followed by 40 cy- 
cles of 95°C for 15 s, 58°C for 30 s and 72°C for 30 s. In 
order to control the specificity of amplification, the reac- 
tion was concluded with a melting curve analysis ramp- 
ing from 55°C to 99°C in steps of 1°C every 5 s. PCR 
efficiencies and cycle thresholds were obtained using 
LinRegPCR 12 [67] and differential expression ratios 
were calculated by the Cx difference formula [68]. Tubu- 
lin beta chain (CC1G_04743) was used as a house keep- 
ing normalizer. In addition, water or 1 ng/uL RNA were 
included as negative control reactions. To further valid- 
ate the significance of the RNA-seq-derived differential 
expression analysis, the constitutive expression of an 
array of housekeeping loci commonly used in qRT-PCR 
normalization [69-75] was verified in the sequencing 
datasets after library size normalization. 

Bioinformatic analysis 

Fastqc files were used to filter and trim the reads to be 
mapped. Strand-specific reads were mapped to the third 
annotation (September 2010) of the C. cinerea oka- 
yama7#130 genome hosted at the Broad Institute (C. 
cinerea Sequencing Project, Broad Institute of Harvard 
and MIT (http://www.broadinstitute.org/)). SOLiD mapped 
reads were counted using Cufflinks [76] . All the sequences 
were deposited in the ArrayExpress database (http://www. 
ebi.ac.uk/arrayexpress/) under the accession number [E- 
MTAB-1968]. To determine the percentage of loci 
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showing baseline expression, 5 reads/locus were taken as 
a minimal threshold [77] . Library size normalization (scal- 
ing method), fold change calculation per locus and fisher 
exact test comparing mycelia and primordia libraries were 
performed using the edgeR package [78] using sense-read 
counts. Fisher exact test p-value < 0.05 and fold change 
> 8 were the criteria established to classify a locus as dif- 
ferentially expressed in one of the two compared samples. 
Functional annotation of loci found to be developmentally 
regulated was explored by PSI-BLAST [79]. Enrichment 
of annotation terms in differentially expressed loci was 
visualized in tag clouds constructed by Wordle (© IBM 
Corporation) [32]; common tags such as "protein", "hypo- 
thetical", "domain", "containing", "family", "putative", 
"similar", "probable", "related" and "subunit" were ex- 
cluded. A complementary functional annotation clustering 
was performed with the Database for Annotation, 
Visualization and Integrated Discovery (DAVID) 6.7 [33] 
using default annotation categories and the C. cinerea 
okayama7#130 translated genome as background. SignalP 
4.1 [80] and TMHMM v. 2.0 [81] were used to predict the 
presence of secretion signal and transmembrane helices in 
developmentally regulated loci. Presence of nuclear 
localization signals (NLSs) in members of the Velvet pro- 
tein family was assessed with NLStradamus [82]. To 
analyze the similarity in the gene expression programs as- 
sociated to early stages in fruiting body development 
across different basidiomycetous fungi, orthologous genes 
to differentially expressed loci in C. cinerea SIP were iden- 
tified by PSI-BLAST (best hit showing an E-value < 0.005) 
in the genomes of Schizophyllum commune 4-39/4-40 and 
Laccaria bicolor S238N-H82. Expression data were re- 
trieved from the original studies where the gene expres- 
sion of vegetative mycelium and SIP or YFB from S. 
commune [43] and L. bicolor [44], respectively, was mea- 
sured. First, log2 expression ratios were calculated for 
stage 1 primordia against monokaryon mycelium from S. 
commune using tags per million per locus. Similarly, log2 
ratios were calculated for "young fruiting bodies" com- 
pared to "free living mycelium" (the latter corresponding 
to the average of two replicates) of L. bicolor from 
quantile-normalized Robust multichip average values. A 
centroid-linkage hierarchical clustering analysis of the 
computed log2 ratios was performed with Cluster 3.0 
and heat maps were generated and visualized using 
TreeView 3 [83]. 

A more detailed orthology and expression analysis was 
carried out for loci known to encode proteins interacting 
with Velvet domain-containing proteins involved in sex- 
ual development and secondary metabolite production 
in Aspergillus [42]. The sequences of RosA, NosA, StuA, 
NsdD, PpoA, VelB, LaeA, KapA FphA and CryA from A. 
nidulans were used to retrieve the corresponding ortho- 
logs from the protein-coding genomes of C. cinerea 



Okayama 7, S. commune 4-39/4-40 and L. bicolor 
S238N-H82. Orthologs were defined as the best hits in a 
PSI-BLAST showing E-values < 0.005. Progressive mul- 
tiple sequence alignments for the orthologs present in A. 
clavatus, A. fumigatus, A. oryzae and A. nidulans, as 
well as in the three basidiomycetous species considered, 
were constructed in CLUSTAL W [84]. Finally, Velvet 
domain-containing proteins in C. cinerea, S. commune 
and L. bicolor were identified in Pfam and SMART dia- 
grams were created [85]. 

LC-MS/MS based shotgun proteomics 

Twenty mg lyophilized VM from a single plate or a SIP 
pool were lysed in a FastPrep FP120 during 40 s at 6 m/s 
using 200 mg 0.5 mm glass beads and 600 uL PBS sup- 
plemented with 1 mM phenylmethanesulfonylfiuoride 
(PMSF) and lx Complete Protease Inhibitor Cocktail 
(Roche). The lysates were centrifuged for 15 min at 16000 
xg and 4°C, and the supernatants (soluble protein ex- 
tracts) recovered. Soluble protein extracts from VM and 
SIP were further processed by filter-aided sample prepar- 
ation (FASP) method as previously described [86] . In brief, 
protein extracts were loaded onto an Amicon column 
equipped with a 10 kDa MWCO membrane (Millipore), 
reduced with 55 mM dithiothreitol at 37°C for 1 h and 
alkylated with 65 mM iodoacetamide in the dark at 37°C 
for 1 h. Reduced and alkylated extracts were digested with 
sequencing grade porcine trypsin (Roche) for 18 h at 37°C 
in 25 mM ammonium bicarbonate, pH 8.5. Digested pep- 
tide mixtures were collected by centrifugation and dried in 
a Savant SpeedVac (Thermo Scientific). All samples were 
de-salted by C18 ZipTip before mass spectrometry analysis. 

Samples were analyzed on a LTQ-Orbitrap Velos mass 
spectrometer (Thermo Fischer Scientific) coupled to an 
Eksigent-Nano-HPLC system (Eksigent Technologies). 
Peptides were suspended in 2.5% acetonitrile and 0.1% 
formic acid, loaded on a self-made tip column (75 um x 
80 mm) packed with reverse phase C18 material (AQ, 
3 um 200 A, Bischoff GmbH) and eluted with 250 nL/ 
min flow rate in a gradient from 3% to 50% of B in 
90 min, 97% B in 10 min. One scan cycle comprised a 
full scan MS survey spectrum, followed by up to 20 se- 
quential CID MS/MS on the most intense signals above 
a threshold of 1500. Full-scan MS spectra (300-2000 m/z) 
were acquired in the FT-Orbitrap at a resolution of 60000 
at 400 m/z, while CID MS/MS spectra were recorded in 
the linear ion trap. CID was performed with a target value 
of 1E4 in the linear trap, collision energy at 35 V, Q value 
at 0.25 and activation time at 30 ms. AGC target values 
were 5E5 for full FTMS scans and 1E4 for ion trap MSn 
scans. For all experiments, dynamic exclusion was used 
with one repeat count, 15 s repeat duration, and 60 s exclu- 
sion duration. For quantitation, each sample was measured 
in technical triplicates using the same parameters. 



Plaza et al. BMC Genomics 2014, 15:492 
http://www.biomedcentral.com/1471-2164/15/492 



Page 14 of 17 



All MS data from VM and SIP were converted to a 
peak list and searched against the C. cinerea oka- 
yama7#130 proteome (predicted transcript translation) 
hosted at the Broad Institute (C. cinerea Sequencing 
Project, Broad Institute of Harvard and MIT (http:// 
www.broadinstitute.org/)) using the Mascot search en- 
gine (version 2.3) considering variable modifications: 
carbamidomethylation on cysteine and oxidation on me- 
thionine. The tolerance of mass accuracy of MS and 
MS/MS was 8 ppm and 0.6 Da. The false discovery rate 
of proteome dataset was 1% and the score of each pro- 
tein exceeded 30. 



Cloning and recombinant expression of CC1G_11805- 
encoding cDNA 

cDNA encoding CC1G_11805 was amplified from C 
cinerea SIP-derived cDNA using the forward and reverse 
primers 5 -CCAGCTTAAAGGAGTCACAAGG-3 ' and 
5'-AACGTTCAACGCCCAGCCAC-3', respectively, with 
a Pfu DNA polymerase. 3' adenines were added to the 
PCR fragment before ligation to the DNA amplification 
vector pGEM-T Easy (Promega). This pGEM-CClG_ 
11805 construct was used as a template to add Ndel and 
BamHI restriction sites using the primers 5-GGGGG 
GCATATGTCTCAAGCAGGGATCACAC-3 ' and 5 -GG 
GGGGGGATCCTCAGATACGCCCGATGACTTC-3', 
respectively. The resulting product was digested with Ndel 
and BamHI in 2x Tango buffer (Thermo Scientific), li- 
gated to a pre-digested pET24 expression vector (EMD 
Millipore), and used to transform chemocompetent E. coli 
BL21. Transformed colonies were selected on LB plates 
containing 50 ug/mL Kanamycin. CC1G_11805 recom- 
binant expression and solubility were assessed as follows: 
pET24-CClG_l 1805-containing E.coli BL21 was grown in 
LB broth supplemented with 50 ug/mL Kanamycin up to 
OD 60 o: 0.5; thereafter, the cultures were divided in two 
separated flasks and protein expression was induced for 
18 h at 24°C by adding 1 mM Isopropyl [3-D-l-thio- 
galactopyranoside (IPTG, final concentration) to only one 
of the two subcultures. Cells were harvested, suspended in 
1 mL PBS containing with 1 mM phenylmethanesulfonyl- 
fluoride (PMSF) and lysed in a single FastPrep FP120 
homogenization step of 35 s at 6 m/s in the presence of 
1 g 0.1 mm glass beads. The whole cell protein extract 
(WCE) was centrifuged for 5 min at 5000 xg and 4°C (low 
speed centrifugation, LS) to remove cell debris and the 
supernatant transferred to a fresh tube to be centrifuged 
again for 30 min at 14000 xg and 4°C (high speed centrifu- 
gation, HS). Samples were collected from the whole cell 
extracts (WCE and WCE + IPTG) and the two centrifuga- 
tion supernatants (LS and HS), mixed with Laemmli buf- 
fer and denatured at 95°C for 5 min to be loaded on a 12% 
SDS-PAGE. 



Caenorhabditis elegans toxicity assay 

In order to evaluate the developmental toxicity of 
CC1G_11805 to nematodes, a C. elegans toxicity assay 
was performed as previously described [40]. In brief, C. 
elegans N2 and C. elegans pmk-1 (km25) (kindly pro- 
vided by M. O. Hengartner) were grown on NGM plates 
(50 mM NaCl, 2.5 g/L bacteriological peptone, 13 mM 
cholesterol and 1.7% agar) pre-seeded with E. coli OP50. 
Eggs were obtained by bleaching gravid hermaphrodites 
in 15 mL conical bottom tubes using a solution contain- 
ing 0.5 N NaOH and 1% NaCIO for 10 min and washing 
twice in 10 mL distilled deionized water (ddH 2 0). Clean 
eggs were transferred to a 1.5% agar plate and hatched 
for 18 h at 20°C. LI larvae were collected in PBS, 
counted and adjusted to 1.5 larvae/uL. Approximately 
30 LI larvae/replicate (4 replicates/treatment) were 
mixed in 200 uL PBS with OD 600 : 2 IPTG-pre-induced 
E. coli BL21 expressing CC1G_11805 in flat bottom 96 
well-plates. As a negative control, IPTG-pre-induced E. 
coli BL21 cells transformed with empty vector pET24 
were used. Worms were incubated for 48 h at 20°C and 
the number of individuals reaching each developmental 
stage or dying was counted. A Mann Whitney test be- 
tween the empty vector pET24 control and CC1G_ 
11805 was run to test the statistical significance of the 
results observed. 

Aedes aegypti toxicity assay 

CC1G_11805 toxicity towards A. aegypti Rockefeller 
(kindly provided by W. Rudin and P. Miiller) was tested 
as previously described [40]. In brief, 600-800 eggs were 
hatched in a glass petri dish containing 200 mL deion- 
ized water and 30 mg ground fish food for 20 h at 28°C 
in the dark. Larvae were transferred to fresh 800 mL de- 
ionized water supplemented with 50 mg ground fish 
food and incubated for 10 h at 28°C in order to obtain 
synchronized L2 larvae. Ten L2 larvae/replicate (4 repli- 
cates/treatment) were starved in 100 mL fresh deionized 
water in 100 mL Schott flasks for 6 h at 28°C before 
adding 1 mL OD 600 : 20 (Final OD 600 : 0.2) IPTG-pre- 
induced E. coli BL21 (empty vector pET24, pET24- 
CGL2 [18] or P ET24-CC1G_11805). Larvae were 
incubated for 96 h at 28°C in the dark and the surviving 
individuals were counted. Percentages of surviving lar- 
vae/treatment were calculated. A Dunn's multiple com- 
parison test between the empty vector pET24 control 
and the different treatments was run to test the statis- 
tical significance of the results observed. 

Availability of supporting data 

The data sets supporting the results of this article are 
available in the ArrayExpress repository under the ac- 
cession number E-MTAB-1968 (http://www.ebi.ac.uk/ 
arrayexpress/). 
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Additional file 1: Table SI. Genome-wide differential expression 
analysis of C cinerea SI P and VM. 

Additional file 2: Table S2. Genes up-regulated in SIP and VM. 

Additional file 3: Figure SI. RNA-seq data validation by qRT-PCR. The 
expression of four selected genes was validated by qRT-PCR. A comparable 
relative expression pattern was found for all the genes evaluated in SIP 
and VM with both techniques. RNA-seq data corresponds to the mean 
log2(Sl P/VM) of two biological replicates. qRT-PCR data show the mean 
log2(Sl P/VM) of three technical replicates from a single biological replicate 
of SIP and VM. Bars correspond to standard deviations. Dashed lines 
indicate the differential gene expression thresholds selected for this study 
(log2(S!P/VM) = +/-3). 

Additional file 4: Table S3. RNA-seq expression of reference 
housekeeping loci. 

Additional file 5: Table S4. DAVID functional annotation enrichment of 
genes up-regulated in SIP and VM. 

Additional file 6: Figure S2. Scheme of fruiting body development in 
C. cinerea. Gene mutations preventing fruiting body development at 
different stages in C cinerea are shown in blue. White and black sections 
indicate light and dark periods, respectively, corresponding to 12 h each. 
Numbers indicate the measured fold (SI P/VM). A mutation in the 
SI P-specific locus cfsl stops sexual development at the initials stage. 
Figure adapted from Kues U, 2000. 

Additional file 7: Figure S3. Differential transcription comparison 
between genes coding for Velvet domain-containing proteins from three 
different basidiomycetes. Amino acid sequences corresponding to the full 
set of Velvet domain-containing proteins encoded in the genomes of C. 
cinerea, L bicolor and S. commune were retrieved from Pfam (PF1 1754). 
SMART domain architecture diagrams are shown in (A). (B) Loci expression 
in SI P (C. cinerea and S. commune) or YFB (£.. bicolor) relative to VM for all 
the velvet domain-containing proteins in the three basidiomycetes 
compared. N/E: No expression detected. 

Additional file 8: Table S5. Orthology and expression of C. cinerea, 
S. commune y L. bicolor genes. 

Additional file 9: Multiple sequence alignment between Velvet- 
interacting proteins of C. cinerea, L. bicolor, S. commune, A. clavatus, 
A. fumigatus, A. oryzae, A. niger and A. nidulans. 

Additional file 10: Figure S4. Differential transcription comparison 
between genes coding for Velvet-associated proteins from three different 
basidiomycetes. C. cinerea, L bicolor and 5. commune loci homologous to 
those encoding VelvetA-associated proteins in the ascomycete Aspergillus 
clavatus (shown on top) were identified by PSI-BLAST (Best hit showing 
an E-value < 0.005). Differential expression in SI P or YFB relative to VM is 
shown for the orthologs identified. log2(SlP or YFB/VM) > 0 indicates 
increased expression in SIP (C. cinerea or 5. commune) or YFB (L bicolor). 
On the contrary, a log2(Sl P or YFB/VM) < 0 represents a decreased 
expression in SIP or YFB. Blue and white locus IDs differentiate neighboring 
groups of orthologs in the chart. Expression of velB and kapA is conserved 
during sexual development among basidiomycetes. N/E: No expression 
detected. 

Additional file 11: Figure S5. Differential transcription comparison 
between genes coding for transcription factors involved in sexual 
development in 5. commune from three different basidiomycetes. 
Orthologs of transcription factors involved in sexual development in 
S. commune were identified by PSI-BLAST in the genomes of C cinerea 
and L bicolor. log2((SlP or YFB)/VM)<0: Down-regulation in SIP or YFB; 
log2((S!P or YFB)/VM)>0: up-regulation in SIP or YFB. Domain architecture 
for each ortholog in the three species compared is shown in the right hand 
panel. 

Additional file 12: Table S6. LC-MS analysis of C. cinerea SI P and VM. 

Additional file 13: Table S7. Genes up-regulated in SI P and VM 
detected by LC-MS. 

Additional file 14: Table S8. qRT-PCR validation primers. 
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